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^ I In quenched QED we construct a non-perturbative fermion-boson vertex that ensures the 

r-| ■ fermion propagator satisfies the Ward-Takahashi identity, is muhiphcatively renormahz- 

able, agrees with perturbation theory for weak couphngs and has a critical coupling for 

P-i. 

QJ . dynamical mass generation that is strictly gauge independent. This is in marked con- 

trast to the rainbow approximation in which the critical coupling changes by 50% just 
K^ . between the Landau and Feynman gauges. The use of such a vertex should lead to a 

c^ ' more believable study of mass generation. 



1 Introduction 

The standard model is highly successful in collating experimental information on the basic 
forces. Yet, its key parameters, the masses of the quarks and leptons, are theoretically 
undetermined. In the simplest version of the model, these masses are specified by the 
couplings of the Higgs boson, couplings that are in turn undetermined. However, it could 
be that it is the dynamics of the fundamental gauge theories themselves that generate the 
masses of all the matter fields. To explore this possibility, the favourite starting point is 
to consider quenched QED [1-9] as the simplest example of a gauge theory and study the 
behaviour of the fermion propagator. Then in the rainbow approximation, it is well-known 
that the fermion field can have a dynamically generated mass if the interaction is strong 
enough, i.e. the coupling, a, is larger than some critical value, ac- This critical coupling 
marks a change of phase and so its value should be gauge independent. Unfortunately, 



the rainbow approximation allows a far from gauge invariant treatment [0, |TT|. The 
purpose of this paper is to construct a non-perturbative fermion-boson interaction that 
respects the Ward-Takahashi identity, ensures the fermion propagator is multiplicatively 
renormalizable, agrees with perturbation theory when a -C 1, and possesses a gauge 
independent critical coupling. 

The gauge technique of Salam, Delbourgo and collaborators [|T^ was developed to 
solve essentially such constraints. However, despite formal results on the first two of 



these [|13|, |lj], their expression in terms of the spectral representation for the fermion 
propagator has proved difficult in practical calculations of the fermion propagator, for 
example, 0. Consequently, we develop an explicit construction procedure amenable to 
straightforward computation. 

We start with the Schwinger-Dyson equation for the fermion propagator. This non- 



linear integral equation encodes all we can know about the fermion propagator. To be able 
to consider this equation alone of the infinite set of Schwinger-Dyson equations - one for 
each Green's function - we must make an ansatz for the full fermion-boson vertex. Quite 
generally, this vertex can be regarded as the sum of two components : the longitudinal 
and transverse parts. The well-known Ward-Takahashi identity constrains the longitu- 
dinal part. How to fulfill this constraint in a manner free of kinematic singularities has 



been solved some time ago by Ball and Chiu |T5|. That multiplicative renormalizability 



constrains the transverse vertex has also been known for some time ||16| , |1J|. However, it 
is more recently that Curtis and one of the present authors |l^ explicitly constructed a 
simple form (perhaps the simplest possible form) to ensure the multiplicative renormal- 
izability of the fermion propagator. This ansatz is called the CP vertex. 

Subsequent study has shown that with this vertex the fermion propagator still has 



the possibility of a chiral symmetry breaking phase [|I8| . Moreover, in dramatic contrast 
to the rainbow approximation, the critical coupling required is only very weakly gauge 
dependent in the neighbourhood of the Landau gauge. However weak this variation, 
any gauge dependence shows that the CP vertex cannot be the exact choice. Here, we 
determine the constraints on the full fermion-boson vertex that ensures gauge covariance 
for the fermion propagator and exact gauge independence for the critical coupling. The 
resulting vertex involves two unknown functions Wi and W2-, which each satisfy a sum 
rule and a constraint on their derivatives. Any choice of these fulfills our fundamental 
constraints as long as it correctly matches onto perturbation theory. This construction 



builds on the CP vertex, extending the work of Dong et al. [T9|. Though the discussion 



in Sect. 2 of how to ensure the gauge covariance of the wavef unction renormalization of 



the fermion propagator is very close to that of Dong et al. ||19[, to make the extension to 



the gauge independence of the critical coupling clear, we have given all the details of our 



formulation making our construction in Sect. 3. self-contained. 

In general, only the position of the pole in a propagator has to be gauge independent. 
At that value of the momentum, when p^ = im? in Minkowski space, (or equivalently at 
p^ = —m? in the Euclidean space in which we work) the fermion mass function has to 
be independent of the gauge. Atkinson and Fry |2^ proved this independence follows 
from the Ward-Takahashi identities. However, at the critical coupling for dynamical 
mass generation, multiplicative renormalizability imposes such a simple form on the mass 
function that this whole function becomes gauge independent. This is embodied in our 
construction. 



Our results have to be compared with earlier work. For example, Rembiesa [£ll and 
Haeri [^, using the previously mentioned gauge technique, construct fermion-boson ver- 
tices that make the fermion propagator itself gauge independent. This is, of course, at 
variance with its behaviour in perturbation theory and consequently with the renormaliza- 



tion group in the weak coupling limit. Rembiesa [21| then went on to find that the critical 



coupling for mass generation with such a vertex is strongly gauge dependent, being given 
by oic = 7r/(3 + ^). In complete contrast, Kondo [Q finds a gauge independent coupling 
as here, but at the expense of using a vertex that has singularities. The construction 
presented here aims to overcome these deficiencies. 

2 The Fermion Equation 

The Schwinger-Dyson equation for the fermion propagator, Sf{p), in QED with a bare 
coupling, e, is displayed in Fig 1, and is given by: 

^S^\p) = ^Sf\p)-e'J^^YSF{k)T^{k,p)A,^{q) , (1) 



where q = k — p and Sf{p) can be expressed in terms of two Lorentz scalar functions, 
F{p^) the wavefunction renormahzation and M.{p'^) the mass function, so that 

^^^P> - i>-M{p-) ■ 
The bare propagator Sp{k) = 1/(|^ — mo), where ttlq is the constant (bare) mass. In 
quenched QED, the photon propagator is unrenormahzed and so is given by its bare 
form : 

where the transverse part, A^^,(g), is defined by this equation and ^ is the standard 
covariant gauge parameter. r'^(/c,p) is the full fermion-boson vertex that must satisfy the 
Ward-Takahashi identity 

g%(A;,p) = ^/(fc)-5/(p) . (2) 

We can simplify Eq. (1) by making use of the Ward-Takahashi identity, Eq. (2) : 

Sp\p) = S'p-\p) + ^e' j^,^,l'S,{k)T'^{k,p)KM) 



d^k 4 
d^k 4 



- i^'( I -^^ j SAk)S-Ap) . (3) 

The third term on the right vanishes Q, as it is an odd integral, and we are left with 

To solve this equation we must make an ansatz for the full vertex, T^{k,p). Our 

aim is to construct a vertex that automatically embodies as much of the physics of the 

^This was not noted in Ref. fl^ as pointed out in [p9J, who better remembered Rcf. 
20 of PI than the authors ! 
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interaction as possible. Following Ball and Cliiu ||T5|, we first write the vertex as a sum 
of longitudinal and transverse components : 

T'^{k,p) = TUk,p) + T>i.{k,p) . (5) 

To satisfy Eq. (2) in a manner free of kinematic singularities, which in turn ensures the 
Ward identity is fulfilled, we have (following Ball and Chiu) : 

r^(A;,p) = a(A;2,p2)^/. + h{k\p'){}t+ ^){k + pY - c^k^p^Yk + pf (6) 

where, 

a(fc^p2) = 

Kk\p') = ^ (ttt^-ttt^t llT— ^ . (7) 

and 

g^r^(fc,p) = o , r^(p,p) = o . (8) 

Ball and Chiu wrote down a set of 8 basis vectors Tf^[k,p) for the transverse part [|T5|] , 
that ensures these conditions, Eq. (8), are fulfilled : 

ri;.{k,p) = J2n{k^p\q')Tt{k,p) (9) 

i 

provided that in the limit k —^ p, rj(p^,p^,0) are finite. Our aim is to determine the 

full vertex by requiring the multiplicative renormalizability of the fermion propagator 
and the gauge independence of the chiral symmetry breaking phase transition. Since the 
longitudinal part of this vertex is specified, Eq. (6), this amounts to determining the 
transverse part and hence the Tj of Eq. (9). Of the eight basis vectors, Tj^, four have even 
numbers of gamma matrices and four have odd numbers. 
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It is here that we make three simphfying assumptions. Firstly, we demand that a 
chirally-symmetric solution should be possible when the bare mass is zero, just as in 
perturbation theory. This is most easily accomplished if only those transverse vectors 
with odd numbers of gamma matrices contribute to Tj<{k,p). Then the sum in Eq. (9) 
involves just z = 2, 3, 6 and 8. The corresponding vectors are : 

T^[k, p) = {p^k ■ q - k^p ■ q){f^+ fi) 
n{k,p) = q'Y-q'^ 

m^v) = r{e-p')-{k+pnf^-^) , 

Ttik,p) = -Yp'k^cr.p+p'' fi-k^' ^ , (10) 

where a^^ = ^ [7/.,7i/] • 

The second assumption is that the functions, r,, multiplying the transverse vectors, 
Eq. (9), only depend on k"^ and p^, but not g^. This allows the angular integrations in 
Eqs. (1,4) to be performed. Thirdly, we assume that, in the Landau gauge, the transverse 
component of the vertex vanishes. This is motivated by its large momentum behaviour 
in perturbation theory. There for k"^ ^ p^, the leading logarithmic behaviour is [|17| : 



r^fc,p ^--^ In- 



' k^ 



;ii) 



where as usual a = e^/in. 

The fermion propagator is determined by the two functions F{p'^) and M.{p^). We 
can project out equations for these by taking the trace of Eq. (4) having multiplied by ^ 
and 1 in turn. 



On Wick rotating to Euclidean space, 



F(j92 



" ' 'd'k. ^'*'' 



Att^P^J k^ + M^{k^) g2 



I a{k^p') 



2k -p- — -2A;V + {k' +p)k-p 



2kY + {k"^ + P^)k ■p-—{k'^- P^fk ■ p 



p^ + k • p (/c^ — p'^){k.p — p^) 



+ A^(A;2)c(A;2,p2) 

—^Wr^ [p^ik^-k-p) + M{e)M{p^){k-p~p^ 

+ T2ik^,p^)[ik''+p^)iky-ik-py) 

+ TG{e,p^)[{e-p^)3k-p 
+ Ts{k\p^)\-2ep^ + 2{k-p)^ 



and 



M{p^ 



mo 



" 'd'k ^<*'' 






- h{k\p')M{k 



{k + pr--,{k'-pr 



+ c{k\p') 



(A;' + A;-p)-- (A;^ -p^)(A;^ - A; -p) 



M{p'){e - k ■ p) - M{k'){p ■ k - p' 



(12) 



g2F(p2) 
+ r2(A;2,/)A^(A;2)[-2A;V + 2(A;-p)2l + 3T^{k\p')M{e) 



+ r6(A;^p^)7W(A;2) 3(A;2-p2) 



(13) 



We are only interested in solving this equation when the bare mass, tjiq is zero. One 
solution of the mass equation, Eq. (13), is, as anticipated, M.ijP') = 0. We first consider 
the wavefunction renormalization, F{p'^), in this case. 



Carrying out the angular integrations in Euclidean space gives : 



F(p2 



1 + 



a^ f^' dp F{P) 



Arc 


Jp2 k^ F{p^) 


3a 


j-p'de e e+p^ 

Jo p^ p^ k^ — p^ 


IQn 


3a 


rA^ dk^ k^+p" / 



F{p') 
F{k^) \ 



167r 7p2 A;2 f^2 _ ^2 y F{p'^) 

Stt Jo p^ p^ 

F{e)K,{k\p') , 



a 

8^ 



A;2 



(14) 



where 



K^{k\p') = {k'-3p' 



K2{k\p') = {p'-3k' 



Ts{e,p') + Ts{e,p')- 

+3{e-p')T,{e,p') 

'T,{k\p')+T,{k\p')^ 

+ 3{e-p')T,{e,p') 



\{e+p')r,{ey) 



^{e+p')Me,p') 



(15) 



(16) 



The following treatment turns out to be very close to that of Dong et al. |T^], in a more 



suitable form for our extension to dynamical mass generation. It is convenient to define 
the combination r of r2, r^ and rg, 

Then 

K,{e,p') = {k^-3p')T{k\p') + 3{k'~p')T,{k\p') (18) 

K2{k\p') = {p'-3e)T{e,p') + 3{e-p')T,{e,p') , (19) 

which can be re-expressed in terms of functions with definite symmetry properties when 
k ^^ p. Thus 



1,2 2\ 



K,{ey) = hs{ey) + K{k\p 



(20) 



K2{ey) = hs{k'y)-ha{ey) 



(21) 



where hs{k'^,p'^) and ha{k'^,p'^) are symmetric and antisymmetric respectively under the 
interchange of k and p, 



hs{k\p^) = -[e+p')r{ey)+?>{k^-p^)T,{k\p') 

ha{e,p') = 2{e-p^)T{ey) 



(22) 
(23) 



As discussed in [^, |1^, |1^], multiphcative renormahzabihty requires that the solution of 
this integral equation for the wavefunction renormalization, -F(p^), must be of the form, 

F{p')=A{fY (24) 

gauge covariance requires, z/ = a^/Aii . Burden and Roberts |2^ noted 



As shown in fl9 



numerically that the fermion equation with the simple CP-vertex correctly generates this 



behaviour, even though the authors of Ref. ||2^, |T^ found v = 2a^/(8Tc + a^) as a result 
of not imposing translational invariance on their loop integrations, as discussed earlier. 

This simple power behaviour is generated by the 1 and the first integral on the right 
hand side of Eq. (12). This requires, as noted in Refs. |26, 19|, a cancellation among 
the remaining integrals. Thus multiplicative renormahzabihty imposes the following con- 
straint : 

3 fP'de e e+p" ( F{k^) \ 



2 Jo p'^ p^ k'^ — p^ 



F(p2 



3 r^' de e+p^ / F{e) \ 
+ r^%%F{e)(K{ey) + K{ey 

Jo p^ p^ ^ 

+ /f ^ F{k') {K{k\p') - Kik^p')) = 



(25) 



where F{p ) = A(p Y and the artificial cut-off. A, can be taken to infinity with impunity. 
The scale invariance of the integrals makes it convenient to introduce the variable x, where 
for < k'^ < p'^, X = k'^/p'^, and for p'^ < k'^ < oo, x = p^ /k"^ [^, [T^. Then 
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2 Jo x-1 ^ ' 



+ / dx x"^^ F {p^) (hs{xp'^,p'^) + ha{xp'^,p'^ 

+ / dxx-''-^F{p^) (hs [p^lx,p^) - ha {f/x,p^y) = (26) 



where 

ri(x) = x{l — x'^) — x~^{l — x~'^) 
ri (1/x) = — ri(x) 

Since this equation must hold true at all p"^, the integrands cannot be functions of p^ but 
only of X. Thus 

F{p^)hs{xp'^,p^) = hi{x), 
F{jP)ha{xp^,p^) = h2{x) 

defines hi,h2. Then Eq. (26) becomes 

- / dx rUx) + / dx x"'^^ (hi(x) + h2{x)) 

2 Jo X — 1 Jo 

+ [ dxx-'"^{hi{l/x)-h2il/x))=0 . (27) 

Jo 

The original symmetry of the r's under the exchange of k'^ and p'^ translates as follows in 



terms of the x- variable 191 : 



h{l/x) = x^hix) 

h2{l/x) = -X''h2{x) 



In the most compact way, Eq. (27) can be written as : 

/ dxWi{x) =0 , (2^ 

Jo 
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where 



Wi{x) = ^ ^^ n{x) + (x^^^ + x-^) ih{x) + h{x)) . 
Thus this function Wi{x) fixes TQ{k'^,p'^) and the combination T{k'^,p'^), so that 



(29) 



2 ^2n 



r(A; ,p 



T,{k\p') 



1 1 1 






4 k"^ -p"^ si{k'^,p'^) 


1 A;2+p2 / 1 1 \ 


1 A;2+p2 _ 
^ 3 A;2 - p2 ^ 


2 (F-p2)2l^^(p) ^(p2)J 


1 1 1 

1 


r /A-2 


/ V^ / J 


6 k^ — p"^ si(A;2,p2^ 



(30) 



2 „2\ 



T(A;',p- 



(31) 



where 



2 „2\ 



si(A; ,p 



|_ f (i.=) + 1_ f (p=) 



It is the first term in Eq. (31) that is essentially the CP vertex in the massless theory. Note 
the automatic appearance of the difference (F(k'^) — F{p^) j, which Curtis et al. [|17 
conjectured was the non-perturbative generalization of the leading logarithm behaviour 
in lowest order perturbation theory, Eq. (11). Indeed, agreement with this behaviour is 
naturally achieved if Wi — i> in this limit. 

The vertex can only have singularities for good dynamical reasons. It cannot have 
kinematic singularities. A sufficient condition for this is to assume that each of the Tj 



(i=l,8) is free of kinematic singularities. Ball and Chiu ||T5| found that with their choice 
of basis vectors T-^ this is indeed true at one loop order in perturbation theory in the 



Feynman gauge and Kizilersii pSf has shown this in any covariant gauge at this order 
too. In the present non-perturbative analysis that this continues to hold with the Ball- 
Chiu basis vectors is a plausible simplifying assumption. Thus 



A;2 



lim {k' -p')Te{k\p') = 



(32) 
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which requires 

Wi{l) + W[{l) = -6u , (33) 



as found by ||T9|. Perturbation theory demands Wi{x) be 0{a). While the form of the 
coefficient function tq is determined by the constrained function Wi{x), it is only the 
combination r of r2, t^, rg that is so specified. By imposing the gauge independence of 
the critical coupling for mass generation, we will be able to separate these functions as 
we now show in Sect. 3. 
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3 The Mass Function 



While for a < ac, there is only one solution Ai{p'^) = 0, as a -^ ftc a second non-zero 
solution becomes possible. This solution bifurcates away from the other solution. Bifurca- 
tion analysis allows us to investigate precisely when this happens. In the neighbourhood 
of the critical coupling, terms quadratic in the mass function can be rigorously neglected. 
Thus the wavefunction renormalization, F{p'^), is that of the massless theory, Sect. (2), 
and the equation for the mass function, Al(p^), Eq. (13) with mo = 0, linearizes : 



7W(p2 
7^ 



^ 9 » » o r o 



r B -^(^^ 



3a rp' dk^ 
A-K Jo p^ 



+ — 



3a /-A' dk"' 



F{e) 



F(p2 



Att Jp2 A;2 



M{k' 



2.nk'] 



+ 



k^ 



F(p2) ^ 2{k- 



.2 :;^-^(^') 1- 



. p'^ 



F{e) 



F{k^ 



_^(>,(,VA.(/)^ 



3a /-p' dk"^ . w,9n ^/,9 
An Jo p^ 



e 



Ue-3p'Mk\p') 



+p\,(e,p') + {e-p')T,{e,p') 






o 



+ er,{k\p') + {k'-p')T,{ey) 



(34) 



If this equation is to be multiplicatively renormalizable with a gauge independent bifur- 
cation, then this imposes further constraints on the transverse vertex, Tj {i = 2, 3, 6). We 
first work in the Landau gauge, where we assume the transverse vertex vanishes. This is 
motivated by the perturbative result of Eq. (11). 
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Then we have simply : 



^ ' Arc Jo p^ 

3a r^' de 



M{k' 



MiP 



2(P-p2) 
2(A;2-p2) 



M{k^)-M{p^ 



Mik^)-M{p^ 



(35) 



This equation has the multiphcatively renormahzable solution, 



Mik^) = B(k 



(36) 



where Eq. (35) requires, 



Stt ^ 3 1 

= 1 H 1 IT cot ITS = J [S) . 

3a s 1 — s 



(37) 



There are two roots for s between and 1. Bifurcation occurs when the two roots for s 



merge at s = Sc, specified by f'{sc) = 0. This point defines the critical coupling |J, ^, ^ 



ac = 87r/3/(sc). Numerically, ac = 0.933667 and Sc = 0.470966. A little away from this 
critical point the exponent s in Eq. (36) is given by 



Sr. ± 



2/(s. 



\f"{ 



a 
ar 



(3J 



It is only at the bifurcation point that the simple behaviour of Eq. (36) holds at all 
momenta. There, only when the mass is still effectively zero is there just one scale. A, 
for the momentum dependence of Ai{k'^). Multiplicative renormalizability then forces a 
simple power behaviour. Such a multiphcatively renormahzable mass function must exist 
in all gauges. Consequently, the exponent, Sc, must be gauge independent. Moreover, 
dynamical mass generation marks a physical phase change and so the critical coupling, 
ac, must also be gauge independent. Thus the critical values, ac, Sc, found in the Landau 
gauge must hold in all gauges. This is achieved as follows. We recall Eqs. (14,25) : 



F(p2) 



1 + 



at 



e f^' de F{k^ 



Att Jp'i A;2 _p(p2) 



(39) 
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Multiplying this equation by A^(p^) and subtracting it from Eq. (34), we obtain : 



Mip^ 



^ r'4Mie 



F{k^] 
4tt Jo p"^ " ^ ' F{p'^) 



3a fP^ dk^ 
An Jo p^ 



3a r^' de 
+ — 



M{k^ 



p" 



2(fc2-p2 

A;2 / 



M{k'') 1 



2(A;2-p2 



F{k'') ' 



47r Jp2 k"^ 



A.(*-^<*^ 



+ 



F(p2) 2(P-p 



F(p2 

"' M(e) (i - ^(*^ 



F(p2 



p^ 



3a fp' dk"^ 



Att Jo p^ 



■Mik')F{k') 



2(A;2-p2) 



[M{e)-M{p^)g^^ 



2 „2\ 



6 



ik'-3p')T,ik\p 



+ p^T,{k\p') + [k'-p')T,{k\p') 



3a /-A^ ^^2 



47r J»2 A;2 



M{e)F{e 



6 



(p'-3A;^)r2(A;^p^) 



+ A;V3(A:^p^) + (A;^-p^)r6(A;^p^) 



(40) 



In order for the above equation to reduce to Eq. (35), it must be true that : 



^ rp' dk"^ 
3 Jo p"^ 



Mik' 



F{e) 

F{p') 


Jo 2(A;2 -p2) \^ 
Jo p^ 


F{p') J 

Fie) \ 

F{P') ) 

^{e-3p')r,{e,p') 




+ p^r,{k^p^) + {k^-p')Mk\p') 




+i t^o^'^n^') 


^(/-3A;^)r,(A;^/) 




+ k'^Ts{k\p' 


'') + {e-p')r,{e,p') 



(41) 



at all momentum p and in all gauges C,- 
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This equation can be written as follows : 



Jo V 



F{k^ 



ik' 3-«(*^ 



2(A;2-p2) 



A' - 3M{k^) 



1- 



F(Jf] 



p^ 2(A;2-p2 

Jo p'^ 
+ J^^^M{k')F{k')K,{k',p'' 



A;2 



(42) 



where K^^ikr^p ) and K/^{k\p ) can, like Ki{k ,p ) and K2{k ,p ), be expressed in terms 
of functions with definite symmetry properties under the interchange of k and p : 



2 „2> 



^s(/i; ,p 



2 ^2n 



fl-alA; ,F 



{e-p'r - AkV] T,{k\p') + '- {e+p') T,{e,p') 

+3(k'-p^)r,{e,p') 



p'-p')[{e + p')T,{k\p')-6r,{e,p') 



(43) 



so that 

Ks{k\p') = gs{k\p')+ga{k\p') 
K.ik^p") = gs{k\p') - gaik^p") . 

Introducing the variable x as before and knowing that M.{k'^) ~ (A;2)~*'= and F{k'^ 
{k"^)", Eq. (42) becomes, 



rp -^C rp^ ^C 



^s,~l + x'^-''-^ 



Jo 2 Jo X — 1 

- / dxx^'^^Fip^) \gs{xp^,P^)+ga{xp\p^ 

Jo L 

dxx'^-''^^F{p'^)[gs{pyx,p'^)-ga{pyx,p'^]) =0 . 



,2 „2^ 



(44) 



Once again, this equation must hold true for all p^, and so the integrands cannot be 
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functions of p^ but solely oix. Thus, we conveniently define, 
Then we have 



U dxx''-'^ + - / 

Jo 2 Jo X — I 



x-'^ - x^-'^ - x'^-^ + x'^-''-^ 



f\xx''-'^[gi{x)+g2{x)]- C dxx'^-'-^[gi{l/x)-g2{l/x)] = Q . (45) 
Jo JO 



The symmetry of the vertex |19| under k ^^ p means that 



gi{l/x) = x^gii^x) 
g2 (l/x) = -x''g2{x) 



In contrast to our discussion in Sect. 2 when the equations for the wavefunction renor- 
malization, -F(p^), apply for all values of the coupling, Eqs. (44,45) only hold when 



a = ar 



Eq. (45) can be written in a compact way as 



1 dx 



Woix) = 



Wx 



(46) 



where 



Wo(x) 



^j.y-s.+^ + 



3 r2{x) 



1 x-\ 



-x" 



\g\{x) +^2(a;)] 



X 



-i/ + Sc 



-5 \g^{\lx)-g2{\lx)\ 



(47) 



with 



r2\x) 



x^' 



[i-xn-x 



[i-x-n 



(48) 
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which has the property, r2 {l/x) = —r2{x). Conveniently defining the combination, 

s.ie p') = ^ :^^ Fie) + p ^^^ F(p') 



(49) 



we have 



2 J2\ 



9sik ,p 



e 



2s2(A;2,p2) 
3 k"^ +p 



k M{k^)F{k^) p Mip^)F{p^] 



p M{p^)F{p^) k Mik^)F{k^) 
1 fe \ 



A k"^ — p'^ S2{k'^,p^] 



''"' \¥ 



2 S2{k'^,p'^ 



Wo 



p^ 



p" 



(50) 



2 J2\ 



9aik ,p 



e 



2s2(A;2,p2) 
3 1 



k M{e)F{e) p M{p^)F{p^ 



p M{p^)F{p^) k M{k^)F{k^) 



4 S2{k\p^) 

1 1 

2 S2{k'^,p'^] 



'k^ 



Wo 



Wo 



pz 



p" 



(51) 



Solving the last two equations for T2 and Ts in terms of tq and W2-, we obtain 



2 J2\ 



2 J2\ 



T2ik ,P 



2^ q2{k\pn ^r^{k\p 



{k^-p^f S2{k\p^) ^ 
1 1 


(P 
W2 

W2 


-p2) 

/P \ /n2 \ 1 

7)2 / 2 7,2 


e + p^ 1 

(fc2-p2)3 S2(A;2,p2) 



(52) 



where 



2 ^2\ 



?2(/i; ,p' 



A;2 — J92 



A;3 M{k'^)F{k^) p^ M{p^)F{p^ 



p M{p^)F{p'^) k Mik^)F{k^] 



(53) 



where q2{k\p ) is obviously a symmetric function of k and p, and 
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2 J2\ 



Tsik ,p 



k? — p^ 

1 
+ 



r,{k\p') 
1 



^2 _ p2 g^(J^2^p2^ 

1 A;2+p2 1 



2 \p 



.2 



-\q.{k\p') 



6(A;2— p2)2 g^(J^2^p2-^ 

.4 g^2p2 ^ 



"'^ ? + "^^^ fe 



p- 



6 (/c2— p2)3 S2{k'^,p'^) 



''4 



w, ^ 



p 



where 



2 „2\ 



93 (A; ,p 



kp 



(P_p2)2 



(p=^ - 3/t=^ 



M{k^)F{k^\ 
M{p^)F{p^] 



- {e - 3p2 



M{k^)F{k^) 
M{p^)F{p^) 



where qz{k'^-,p^) is antisymmetric in k and p. The relation, Eq. (17), 



r(A;^p^) 



r3(A;^p2) + r8(A:^p^) - i (P +/)r2(P,p2) 



then fixes rs{k'^,p'^). 
rs{k',p') 



-2'^^,T,{e,p')+r{e,p') 

1 1 



A;2 — p2 g^(J^2^p2^ 

1 e + p'^ 1 



^2 



3 (A;2-p2)2 s2(A;2,p2) 
2 A;4 + p4 1 



„2 



Wo 



e 



2 „2\ 



g8(A; ,P' 



pz 



+ Wo 



3 (A;2 — p2)3 g^(J^2^p2 



"'' u^ 



A;2 



where 



(54) 



(55) 



(56) 



(57) 



2 ^2\ 



98 (^ ,P 



1 



(A;2 — p 



2^2 



'^(3^^+/)4S^-?(^^ + 3.-^(^^)^(^^) 



(5^ 



Al(p2)F(p2) A;' ^ ' 7W(p2)F(p2) 

which is clearly antisymmetric in k and p. Imposing the condition that the vertex and its 
components should be free of kinematic singularities means that. 



\im(k'-p')n{ey) = Q ^ = 2,3,8 
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noting that the antisymmetry of rg means Tq{p^,p'^) = 0. Thus 

W2{l) + 2W^{l) = 2aiy-s + l) , (59) 

where s = Sc at the critical point. The transverse vertex has the correct lowest order 
perturbative limit, viz. Fj. = 0{a), provided, 

2/ 2\ . k M{k'^)F{e^ 



Since at large momenta we expect the power behaviour of Eqs. (24,36) even away from 
criticality, Eq. (59) will hold for all values of the coupling, a. In contrast, Eq. (46) is 
only true at the bifurcation point. Its exact form for all a is not known, but Eq. (38) 
might suggest 



'^W,{x) ^eJl-- , (61) 

Wx V ar 



to agree with both the a = and a = ac limits, Eqs. (60,46). These equations determine 
our vertex for any Wi {i = 1, 2) that satisfy the constraints. 



In 1 18 1 a plot is shown of the critical coupling, a^ as a function of the covariant gauge 
parameter, ^, when the CP vertex is used. We refrain from showing the analogous graph 
for the presently constructed vertex, as ac would be boringly gauge independent ! This 
has been achieved for any choice of the functions Wi{x) {i = 1,2), that satisfy Eqs. (28, 
33, 46, 59, 60). A simple example for Wi is 2z/ (1 — 2a;). There are, of course, an infinity of 
such functions. In practice, we expect that Wi should be expressible solely in terms of the 
ratio F{k'^)/F{p'^), while W2 should surely also involve Ai{k'^)/Ai{p'^). However, we have 
not been able to find simple examples that achieve this. The exact form of the full vertex 
would, of course, determine these functions precisely. Thus solving the Schwinger-Dyson 
equation for the three point function would specify the unknowns. However, that has 
not been our aim. Our aim is more limited. It is to construct a vertex that ensures the 
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fermion propagator is gauge covariant, multiplicatively renormalizable and has a gauge 
independent chiral symmetry breaking phase transition. One does not need to know 
the exact form of the full vertex to achieve these properties only the effective vertex 
for the fermion equation, Eq. (1). However, we believe that this effective vertex should 
nevertheless satisfy the appropriate Ward-Takahashi identity and agree with perturbation 
theory at least in the leading logarithmic limit of the weak coupling regime. This is the 
construction we have achieved for any functions Wi{x) {i = 1,2). This effective vertex is 
thus given by Eqs. (5,6,7,9,10,31,51,53,56). 

4 Conclusions 

The non-perturbative behaviour of the fermion propagator is governed by its Schwinger- 
Dyson equation . In quenched QED, the self-consistent solution of this equation is de- 
termined by the fermion-boson interaction. This in turn satisfies a Schwinger-Dyson 
equation that relates it to the full 4-point function and this 4-point function satisfies its 
own Schwinger-Dyson equation While the solution of this infinite set of equations repre- 
sents the whole theory, the complete set is, of course, impossible to solve. Consequently, 
we need a systematic method of truncation that maintains the key features of the theory 
: its gauge invariance and multiplicative renormalizability. The only known truncation 
scheme consistently respecting these properties is perturbation theory. However, the bulk 
of strong interaction phenomena require a non-perturbative approach. Thus, for example, 
massless bare matter fields remain massless to all orders in perturbation theory. However, 
if the interactions are strong enough, a chiral symmetry breaking phase may become a 
possibility. Truncating the nested Schwinger-Dyson equations to just the fermion equa- 
tion by the rainbow approximation, in which the fermion-boson vertex is simply treated 
as bare, this possibility is realized. However, this approximation is highly gauge depen- 
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dent with the critical couphng for this phase transition varying by a factor of two from 



(^ = to 3 [^. The present paper defines a truncation of the fermion Schwinger-Dyson 
equation , which does respect the key properties of the theory. The vertex constructed 
satisfies the Ward-Takahashi identity, ensures the fermion propagator is multiphcatively 
renormahzable, agrees with one loop perturbation theory for large momenta and enforces 
a gauge independent chiral symmetry breaking phase transition. This is a step on the 
way to a meaningful non-perturbative truncation scheme : meaningful in the sense that 
the fundamental aspects of the physics crucially determining the fermion propagator are 
thereby encoded in its Schwinger-Dyson equation. 

Investigation of how, for a given coupling strength, the generated mass compares 
with that found using the rainbow approximation requires the solution of the coupled 
equations for F{p'^) and M.ijP'). Study of the chiral symmetry breaking phase transition, 
using bifurcation analysis, fortunately allows these equations to be uncoupled rigorously. 
The coupled solution is planned. 

The fact that in a (more) realistic version of non-perturbative QED, mass genera- 
tion is possible makes it more, rather than less, likely that such a phase transition has 



been observed in heavy ion collisions [Q. Moreover, it motivates the need for a realistic 
calculation of ti condensates as the source of the electroweak symmetry breaking [pl|] . 
A realistic calculation, of course, requires the unquenching of the theory. This brings 
at once renormalizations of the transverse photon propagator and of the fermion-boson 
coupling. It is this renormalized coupling, which at the corresponding chiral symmetry 
breaking phase transition, is the physical quantity that must be gauge- independent. The 
need to ensure the multiplicative renormalizability of the now coupled photon propaga- 
tor, of the fermion-boson coupling, as well as of the fermion propagator, significantly 
complicates the problem. The fermion-boson vertex (in particular its transverse part) 
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will intimately depend on the photon renormalization function in a non-perturbative way 
not yet understood. 

Thus the complete multiplicative renormalizability of 2 and 3-point functions brings 
not merely greater algebraic but also methodological complexity. The results for quenched 
QED presented here provide the starting point for such an investigation of full QED. The 
solution to this problem will in turn be the starting point for a study of QCD, where 
boson self-interactions, so essential for both asymptotic freedom and confinement, will 
further complicate the analysis whether in covariant or axial gauges. All this is for the 
future. 
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Figure Caption 

Fig. 1. Schwinger-Dyson equation for tlie fermion propagator. 

Tlie straiglit lines represent fermions and tlie wavy line the photon. 
The solid dots indicate full, as opposed to bare, quantities. 
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This figure "figl-l.png" is available in "png" format from: 
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